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Some asteroids eject dust, producing transient, comet-like comae and tails; these are the 
active asteroids. The causes of activity in this newly-identified population are many and varied. 
They include impact ejection and disruption, rotational instabilities, electrostatic repulsion, 
radiation pressure sweeping, dehydration stresses and thermal fracture, in addition to the 
sublimation of asteroidal ice. These processes were either unsuspected or thought to lie beyond 
the realm of observation before the discovery of asteroid activity. Scientific interest in the active 
asteroids lies in their promise to open new avenues into the direct study of asteroid destruction, 
the production of interplanetary debris, the abundance of asteroid ice and the origin of terrestrial 
planet volatiles. 


1. INTRODUCTION 

Small solar system bodies are conventionally labeled as 
either asteroids or comets, based on three distinct proper¬ 
ties; (1) Observationally, small bodies with unbound atmo¬ 
spheres (“comae”) are known as comets, while objects lack¬ 
ing such atmospheres are called asteroids. (2) Dynamically, 
comets and asteroids are broadly distinguished by the use 
of a dynamical parameter, most commonly the Tisserand 
parameter measured with respect to Jupiter (Kresak 1982, 
Kosai 1992). It is defined by 


Tj = — 
a 




aj 


1/2 


cos(*) 


( 1 ) 


where a, e and i are the semimajor axis, eccentricity and 
inclination of the orbit (relative to Jupiter’s orbit) while aj 
= 5.2 AU is the semimajor axis of the orbit of Jupiter. This 
parameter, which is conserved in the circular, restricted 3- 
body problem, provides a measure of the close-approach 
speed to Jupiter. Jupiter itself has Tj = 3. Main belt as¬ 
teroids have a < aj and Tj > 3 while dynamical comets 
(from the Kuiper belt and Oort cloud) have Tj < 3. (3) 
Compositionally, comets are ice-rich small bodies formed 
beyond the snow-line in the protoplanetary disk, while the 
asteroids are ice-free and formed inside it. 

While there is often a reassuring concordance among 
these classification systems, all three are potentially fallible. 
For instance, the ability to detect a low surface brightness 
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coma or tail is, in part, a function of instrumental parame¬ 
ters and observing conditions. The utility of Equation ([T]| is 
limited for objects with Tj very close to 3 because the un¬ 
derlying criterion is based on an idealized representation of 
the Solar system (e.g. Jupiter’s orbit is not a circle, the grav¬ 
ity of other planets is not negligible, and non-gravitational 
forces due to outgassing and photon momentum can be im¬ 
portant). The least useful metric is the composition be¬ 
cause, except in special cases, we have no practical way to 
measure either the composition of a small body and neither 
can we determine its formation location. 

Taken together, the observational and dynamical classifi¬ 
cations suggest a simple two-parameter schematic that use¬ 
fully describes the solar system’s small body populations 
(Figure [^. The four quadrants in Figure Q conveniently 
separate comets (upper left) from asteroids (lower right) and 
distinguish likely defunct comets (lower left) in which there 
is no activity presumably due to the past depletion of near¬ 
surface volatiles (Hartmann et al. 1987), from the newly 
recognized active asteroid class (upper right). 

The working definition employed here is that active as¬ 
teroids are small bodies which 1) have semimajor axis 
a < aj, 2) have Tj > 3.08 and 3) show evidence for mass 
loss, for example in the form of a resolved coma or tail. In 
defining the critical Tisserand parameter separating aster¬ 
oids from comets as Tj = 3.08, rather than 3.0, we avoid 
many ambiguous cases caused by deviations of the real 
solar system from the circular, restricted 3-body approxi¬ 
mation. This definition also excludes Encke-type comets 
(2P/Encke has Tj = 3.02), and the quasi-Hilda comets 
(Tj ~ 2.9 - 3.04). The orbital distribution of the currently 
known active asteroids is shown in Eigure 0. Three ob- 
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Fig. 1.— Empirical classification of small bodies based on the 
Tisserand parameter, Tj, (x-axis) and the presence or absence of 
coma (y-axis). JFC, LPC and HFC are the Jupiter family. Long 
period and Halley family comet sub-types, distinguished by their 
dynamics. From Jewitt (2012). 


jects occupy planet-crossing orbits while the remaining 14 
orbit in the main-belt. 

The active asteroids are remarkable for being an en¬ 
tirely new population located in one of the closest and most 
intensively-studied regions of the solar system. Their activ¬ 
ity is driven by a surprisingly diverse set of physical pro¬ 
cesses. Reported examples include hypervelocity impact 
(e.g. (596) Scheila), rotational instability probably driven 
by YORP torques (e.g. 31 IP, P/2013 R3), thermal disinte¬ 
gration driven by intense solar heating of minerals ((3200) 
Phaethon) and the sublimation of ice ((1) Ceres, 133P, 238P 
and 313P are the strongest examples). Impact and rotational 
disruptions are primary sources of meteorites and larger 
planet impactors. Observations promise to improve our un¬ 
derstanding of the physics, and of the rates, of both aster¬ 
oid destruction processes and meteoroid production. Anal¬ 
ogous disruptions occurring around other stars are likely 
responsible for the creation of debris disks (e.g. Shannon 
and Wu 2011). A common dynamical end-fate of planet¬ 
crossing bodies is to strike the Sun; a better understand¬ 
ing of the role of thermal disintegration will be important 
both in the solar system and in the context of photospheric 
impactors in polluted white dwarf systems (Jura and Xu 
2013). Lastly, the survival of primordial ice in the asteroids 
may offer the opportunity to sample volatiles from a region 
of the protoplanetary disk different from that in which the 
Kuiper belt and Oort cloud comets formed. The outer as¬ 
teroid belt is also a likely source region for the volatile in¬ 
ventory of the Earth (Morbidelli et ai, 2000; O ’Brien et al, 
2006), giving new relevance to the origin of the oceans. 

The active asteroids were called “main-belt comets” by 
Hsieh and Jewitt (2006) while Hsieh et al, (2012a) em¬ 
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Fig. 2.— Distribution of the active asteroids in the semima¬ 
jor axis vs. orbital eccentricity plane. Dynamical asteroids 
are shown as small grey circles, comets as small “H-” sym¬ 
bols, and the active asteroids as black circles, each labeled 
with the object name (c.f Table [T]). Objects plotted above 
the diagonal arcs cross either the orbit of Mars, or of Jupiter. 


ployed the term “disrupted asteroids” to refer to objects that 
exhibit comet-like activity believed to be non-sublimation- 
driven. We use the more general term “active asteroids” 
both because some of the known examples are not in the 
main-belt (c.f. Figure]^, and because this nomenclature im¬ 
plies no supposition about the cause of the activity. 

The active asteroids were reviewed by Bertini (2011) 
and, in more detail, by Jewitt (2012). This chapter updates 
the latter paper with many new observations and examples 
of asteroid activity and adopts a tutorial style in the interests 
of clarity. 

2. CURRENTLY KNOWN ACTIVE ASTEROIDS 

We list the active asteroids in the order of decreasing 
Tisserand parameter in Table Q and discuss them briefly 
here. 

2.1. (3200) Phaethon, Tj = 4.508 

Phaethon is a B-type asteroid with a 0.14 AU perihe¬ 
lion distance, a dynamical lifetime to scattering by the ter¬ 
restrial planets ~100 Myr and a source in the main aster- 
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oid belt that may be related to (2) Pallas (de Leon et ai, 
2010). This 5 km diameter body is the source of the Gem- 
inid meteor stream (e.g. Williams and Wu 1993) and is 
dynamically related to kilometer-sized asteroids 2005 UD 
{Ohtsuka et al. 2006, Jewitt and Hsieh 2006, Kinoshita et 
al. 2007) and 1999 YC (Kasuga and Jewitt 2008). The dis¬ 
persion age of the Geminid stream is t ^10^ yrs {Ohtsuka 
et al. 2006), meaning that Phaethon is active on this or a 
shorter timescale. The timescale for the separation of 2005 
UD and 1999 YC is not known but, presumably, is much 
longer. No gas has been reported in optical spectra {Cham¬ 
berlin et ai, 1996) but near-Sun brightening of Phaethon 
by a factor of two was detected in STEREO spacecraft 
data in 2009 {Jewitt and Li 2010) and 2012 {Li and Je¬ 
witt 2013), while the ejected dust has also been resolved 
{Jewitt et al, 2013, c.f. Eigurej^. However, the sudden 
appearance and position angle of the Phaethon dust tail in¬ 
dicate that the ejected particles are small, with an effective 
radius ~1 pm and a combined mass ~3xl0^ kg {Jewitt 
et al, 2013). This is tiny compared to the Ms ~ 10^^ to 
10^^ kg Geminid stream mass {Hughes and McBride 1989; 
Jenniskens 1994) and suggests that the Geminids are pro¬ 
duced by a different process. In any case, 1 pm parti¬ 
cles are quickly accelerated by solar radiation pressure to 
faster than the solar system escape speed, and cannot con¬ 
tribute to the Geminid stream. Larger particles evidently 
contribute too little to the optical scattering cross-section 
to be discerned in STEREO near-Sun observations taken 
against the bright coronal background. However, particles 
with sizes >10 pm were recently reported in thermal emis¬ 
sion at 25 pm {Arendt 2014), while kilogram-mass Gemi¬ 
nids have been recorded striking the night-side of the Moon 
{Yanagisawa et ai, 2008). Some such bodies might sur¬ 
vive passage through the Earth’s atmosphere {Madiedo et 
al, 2013) and could already be present, but unrecognised, 
in terrestrial meteorite collections. 


Eig. 3.— (3200) Phaethon at perihelion in 2009 and 2012 show¬ 
ing extended emission along the projected Sun-comet line. The in¬ 
sets show the point-spread function of the STEREO camera. Each 
panel shows a region 490”square and is the median of ~30 images 
taken over a 1 day period. Erom Jewitt et al. 2013. 

The Geminid stream mass and dynamical age together 
could imply ejection of debris from Phaethon, if in steady- 
state, at rates 30 < Ms/t < 300 kg s“^. More likely, mass 


loss from Phaethon is highly variable, with dramatic bursts 
interspersed with long periods of quiescence. Continued 
observations of Phaethon, especially at long wavelengths 
sensitive to large particles, are needed. 

2.2. 311P/PANSTARRS (P/2013 P5), Tj = 3.662 

31 IP is an inner belt asteroid (Table[T]l that ejected dust 
episodically over at least nine months in 2013, creating a 
remarkable multi-tail appearance {Jewitt et al, 2013, 2014; 
Hainaut et al, 2014; Moreno et al, 2014, c.f. Eigure |^. 
Interpreted as synchrones (the sky-plane projected positions 
of dust particles of different sizes released simultaneously 
from the nucleus), each tail has a position angle linked to 
the ejection date. The intervals between ejections appear 
random. 
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Eig. 4. — 31 IP observed on two epochs showing the distinctive 
evolving, multiple tail structure. Prom Jewitt et al, (2013). 

The episodic mass loss is unlike that seen in any previ¬ 
ously observed comet. This fact alone argues against ice 
sublimation as the driving agent. An additional considera¬ 
tion is that the orbit of 31 IP lies near the inner edge of the 
asteroid belt, in the vicinity of the Elora family. The Eloras 
have been associated with the LL chondrites {Vernazza et 
ai, 2008), which themselves reflect metamorphism to tem¬ 
peratures ~800°C to 960°C {Keil 2000). It is improbable 
that water ice could survive in such a body. Impact like¬ 
wise offers an untenable explanation for activity that occurs 
episodically over many months. 

The color of 31 IP indicates an S-type classification {Je¬ 
witt et ai, 2013, Hainaut et al, 2014), consistent with its 
inner-belt orbit and with the Eloras. Elora family asteroids 
have a mean visual geometric albedo 0.29±0.09 {Masiero et 
al, 2013). With this assumed albedo, the nucleus of 31 IP 
has a radius r„ < 240±40 m {Jewitt et al, 2013). This 
small size, combined with the inner-belt location, renders 
31 IP susceptible to spin-up by radiation forces. Specifi¬ 
cally, the YORP timescale for 31 IP is < 10® yr, shorter than 
the collisional lifetime. Therefore, it is reasonable to con¬ 
jecture that episodic mass loss from 31 IP results from a ro¬ 
tational instability in which regolith is locally unstable and 
occasionally avalanches off the surface in response to rapid 
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spin (Jewitt et al, 2013) . This “rotational mass shedding” 
qualitatively accounts for the non-steady tail formation and 
the success of synchrone models, which assume ejection 
from the nucleus at zero initial velocity. Unstable mate¬ 
rial on 31 IP would depart at the gravitational escape speed, 
Ve ~ 0.3 m s“^ for a body with p = 3300 kg m“^, the den¬ 
sity of the LL chondrites, or less. Hainaut et al. (2014) 
suggest an alternative model in which dust production re¬ 
sults from friction between two oscillating components of a 
contact binary nucleus. 

2.3. P/2010 A2 (LINEAR), Tj = 3.582 

The object showed a distinctive morphology with a lead¬ 
ing, point-like nucleus about 120 m in diameter (Table |^, 
trailed by an extended tail of dust in which are embed¬ 
ded ribbon-like structures (Jewitt et al. 2010, see Figure]^. 
The position angle of the tail and its variation with time 
are consistent with the action of radiation pressure on mm 
to cm sized dust particles, following impulsive ejection at 
very low speeds (~0.2 m s“^) in 2009 February - March, 
nearly a year before discovery {Jewitt et al., 2010, Snod¬ 
grass et al, 2010). Prediscovery observations were found 
as early as UT 2009 November 22 while detection in the 
first ~6 months after the dust ejection event was impeded 
by the angular proximity of P/2010 A2 to the Sun {Jewitt 
et al, 2011a). Before discovery, a large quantity of fast- 
moving particles are presumed to have left the vicinity of 
the main nucleus. The mass of particles remaining in the tail 
at discovery is estimated to be in the range (6 to 60) x 10" 
kg {Jewitt et al., 2010, Moreno et al., 2010, Snodgrass et 
al, 2010). Observations in 2012 October reveal a surviv¬ 
ing trail of particles up to 20 cm in radius and a differ¬ 
ential power law size distribution index 3.5±0.1 {Jewitt et 
al, 2013b). Kim et al., (2012) found that the colors most 
closely resemble those of an H5 chondrite. The estimated 
total mass is ^5x 10® kg, about 10% of the mass of the nu¬ 
cleus. 



Fig. 5.— Trailing structures in P/2010 A2 observed UT 2010 
January 29. The full image width is ~l' while the inset width is 
~12". From Jewitt et al., (2010). 

The properties of P/2010 A2 appear consistent both with 


an impact origin (with a meter-scale projectile striking at 5 
km s“^) and with rotational instability {Jewitt et al., 2010), 
although an impact origin is often assumed {Snodgrass et 
al., 2010, Flainaut et al., 2012). In the former interpretation, 
the distinctive arms seen in Figure Q are projections of the 
impact cone {Kleyna et al., 2013) while, in the latter, they 
are ribbons of debris shed in the rotational equator {Agarwal 
et al., 2013). 


2.4. (1) Ceres, Tj = 3.309 

Ceres is unique among the known active asteroids in 
that it is the only object for which water sublimation has 
been spectroscopically detected. The 556.936 GHz water 
ground-state line was reported in absorption by Kuppers et 
al., (2014). The line area is time-variable, and can be fit¬ 
ted by a model-dependent water production rate dM jdt ^ 
6 kg s“^. Unlike other active asteroids, however, no dust 
emission has ever been reported for Ceres. 

At the subsolar point on Ceres (heliocentric distance 2.6 
AU) a perfectly absorbing water ice surface would subli¬ 
mate in equilibrium with sunlight at the specific rate = 


5x 10“° kg m“^ s“^ (c.f. Section 3.4 1 . An exposed, subso¬ 
lar ice patch of area {dMjdt)!fs = 0.12 km^, correspond¬ 
ing to a circle of radius 200 m, could supply the measured 
water vapor. A larger area of exposed ice would be needed 
if the ice were more reflective, or if it were located away 
from the subsolar point. One possibility is that ice has been 
recently exposed on Ceres by a surface disturbance, possi¬ 
bly by the formation of a small impact crater or other ge¬ 
ological instability. Alternatively, the water vapor might 
not be produced by sublimation in sunlight, but by subsur¬ 
face heating followed by escape through a vent (c.f. Ence- 
ladus; Porco et al, 2014), since Ceres is potentially large 
enough to maintain significant internal heat {Castillo-Rogez 
and McCord 2010). Additional measurements, presumably 
from the DAWN mission, will be needed to distinguish be¬ 
tween these possibilities. 


2.5. (2201) Oljato, Tj = 3.299 

Magnetometers on the Pioneer Venus spacecraft re¬ 
vealed, in the 1980s, multiple, symmetric interplanetary 
magnetic field enhancements, clumped non-randomly in 
time (Russell et al. 1984). About 25% of these events 
are associated with planet-crossing asteroid (2201) Oljato, 
whose orbit lies interior to Venus’ when near perihelion. 
Russell et al. suggested that the magnetic disturbances re¬ 
sult from deceleration of the solar wind, perhaps caused by 
mass loading from ionized gases released by an unknown 
process from debris distributed along Oljato’s orbit. A mass 
loading rate of only ^5 kg s“^ is reportedly needed. Ob¬ 
servations in 2006-2012 with the Venus Express under a 
similar geometry reveal no events related with Oljato. Lai 
et al., (2014) argue that the field enhancements were due 
to loading of the interplanetary wind by charged nano-scale 
dust in Oljato’s orbit and that the quantity of this dust de¬ 
creased between the Pioneer and Venus Express missions. 
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A spectroscopic search for gas produced by Oljato itself 
proved negative (Chamberlin et al. 1996), with upper limits 
to the CN production rate near 10^^ s“^. With a standard 
H 2 O/CN mixing ratio of 360 (A’Hearn et al, 1995), the 
corresponding limit to the mass production rate in water is 
<1.5 kg s“^. Whatever the cause of the repetitive magnetic 
disturbances, they are not products of an inert asteroid and 
imply mass loss from Oljato. A dynamical simulation in¬ 
dicates that (2201) Oljato has negligible chance of being a 
captured Jupiter family comet (Bottke et al 2002). 

2.6. P/2012 F5 (Gibbs), Tj = 3.228 

P/2012 F5 was observed on 2012 September 18 to ex¬ 
hibit a dust trail extending > 15' in the plane of the sky 
while at a heliocentric distance of i? ^ 3.1 AU. Follow-up 
observations in 2013 showed the object to be largely inac¬ 
tive and set an upper limit on the diameter of the nucleus 
of ~2 km, although residual dust contamination of the nu¬ 
cleus photometry at the time could not be completely ruled 
out (Novakovic et al, 2014). A series of deep-images in 
2014 showed rapid nucleus rotation (period 3.24 hr) and 
revealed four condensations in an orbit-aligned dust trail 
(Figure |^. Dynamical analysis of the object found it to 
be dynamically stable over at least 1 Gyr, and therefore un¬ 
likely to be recently implanted from elsewhere in the solar 
system [Stevenson et al, 2012), although it was also found 
to be a member of an extremely compact asteroid cluster 
determined to be just 1.5 ± 0.1 Myr in age [Novakovic et 
fl/.,2014). 



Fig. 6.— Keck telescope images of P/2012 F5 on UT 2014 Au¬ 
gust 26. The region shown is 1.0x2.5'. The rectangular hox in the 
top panel shows the region enlarged in the bottom panel, where a 
model of the trail has also been subtracted. Letters mark the pri¬ 
mary nucleus, N, and tail condensations, A - D. From Drahus et 
al, (2015). 

The origin of the mass loss is unclear, with evidence con¬ 
sistent both with impact and rotational instability. The dust 
trail was determined via numerical modeling to consist of 
particles ejected in a single impulsive event, consistent with 


impact, roughly nine months prior to the discovery of activ¬ 
ity [Moreno et al, 2012; Stevenson et al, 2012). In contrast, 
the rapid rotation suggests that material may have been lost 
from F5 by rotational instability. Future observations are 
needed to determine the ejection times and fates of the trail 
condensations. 

2.7. 259P/(Garradd) (P/2008 Rl), Tj = 3.216 

259P was observed over a ~45 day interval in 2008 to 
have the appearance of an active comet with a typical flared 
tail while at a heliocentric distance of i? ~ 2 AU (Jewitt 
et al., 2009). The object’s intrinsic brightness decreased 
by a factor of about 2 over the course of those observa¬ 
tions, corresponding to a mass loss rate on the order of 
^ 10“^ kg s“^, assuming mean grain radii of 10 pm and 
bulk densities of p = 1300 kg m“^ (MacLennan & Hsieh, 
2012). Subsequent observations of the inactive nucleus 
found an effective nucleus radius of r = 0.30 ± 0.02 km 
(Table 1^ MacLennan & Hsieh 2012), assuming a red ge¬ 
ometric albedo of 0.05. Spectral observations limited the 
production of the CN radical to Qcn < l -4x 10^^ s“^, cor¬ 
responding to a water production rate <1.5 kg s“^ assum¬ 
ing H 2 O/CN = 360. 259P is located near the 8:3 mean- 
motion resonance with Jupiter and is also affected by the oq 
secular resonance. The dynamical lifetime in this orbit is 
short (20 to 30 Myr) compared to the age of the solar sys¬ 
tem, suggesting that 259P was scattered into its present lo¬ 
cation from elsewhere in the asteroid belt, or possibly even 
from elsewhere in the solar system. 

2.8. (596) ScheUa, Tj = 3.208 

(596) Scheila, a 113 km diameter object with red geo¬ 
metric albedo ~0.04 (Table |^, developed a comet-like ap¬ 
pearance in late 2010 in the form of two prominent dust 
plumes. Over the course of a month, dust in these plumes 
dispersed from the nucleus due to solar radiation pressure, 
apparently without any continued replenishment of parti¬ 
cles from the nucleus [Bodewits et al, 2011; Jewitt et al, 
2011b; Moreno et al, 2011a). Dust modeling by Ishig- 
uro et al (2011a,b) demonstrated that the morphology of 
the observed dust plumes was consistent with the results 
of an impact-driven ejecta cloud consisting of a circularly- 
symmetric ejecta cone (subsequently inverted by radiation 
pressure) and a down-range plume (Figure]^. Bodewits et 
al, (2014) subsequently reported a change in the rotational 
lightcurve which they attributed to the signature of the im¬ 
pact scar. 

Upper limits to the gas production from the nucleus of 
Qoh <10^® s“^ (corresponding to a water production rate 
of < 3 kg s“^; Howell and Lovell, 2011) and Qcn < 
9 X 10^^ (corresponding to QH 2 O ^ 10 kg s“^; Hsieh et 
al, 2012a) were found, though the meaning of these limits 
is unclear given the apparently impulsive nature of the mass 
loss event from Scheila. Jewitt et al (2011) calculated the 
mass of dust in micron-sized grains to be 4 x 10^ kg, while 
more model-dependent attempts to account for larger parti- 
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cles gave total dust masses of ^ 10® — 10^° kg {Bodewits et 
ai, 2011; Ishiguro et al, 2011a; Moreno et ai, 2011). No 
ice was observed in the coma (Yang and Hsieh, 2011). 
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Fig. 7.— Observed image of Scheila on 2010 December 12 (top 
left), simulated image of dust ejecta consisting of an impact cone 
and downrange plume (top right), and diagram of the modeled 
impact scenario (bottom). From Ishiguro et al, (201 lb). 




2.9. 288P/(300163) 2006 VW 139 , Tj = 3.203 

Discovered in 2006 as an inactive asteroid, 288P was 
found to be cometary in 2011 by the Pan-STARRSl sur¬ 
vey telescope (Hsieh et ai, 2012b, c.f. Figure [^. A short 
(~ 10") antisolar dust tail and a longer 60") dust trail 
aligned with the object’s orbit plane were seen in deep 
follow-up images, indicating the simultaneous presence of 
both recent and months-old dust emission, strongly suggest¬ 
ing that the observed activity was due to a long-duration 
emission event, consistent with sublimation. Photometric 
monitoring showed that intrinsic brightness of the near¬ 
nucleus coma remained constant for at least one month, be¬ 
fore then declining by 40% over the next month, again con¬ 
sistent with a long-duration emission event where the coma 
was continually replenished by fast-dissipating small dust 
particles over its period of constant brightness, and then 
faded quickly once the replenishment rate slowed. 

Spectroscopic observations did not detect any gaseous 
sublimation products, setting an upper limit to the CN pro¬ 
duction rate of Qcn < 10^"^ mol s“^ (Hsieh et al, 2012b; 
Licandro et ai, 2013). However, dust modeling found that 
the onset of activity occurred shortly after perihelion, and 
persisted for about 100 days (Licandro et al, 2013). A 
dynamical analysis by Novakovic et al (2012) found that 
288P belongs to a compact cluster of 24 asteroids believed 
to originate from the fragmentation of a ~ 11-km-diameter 
parent body 7.5 ± 0.3 Myr ago. 

2.10. (62412) 2000 SY178, Tj = 3.197 

(62412) 2000 SY178 has diameter 7.8±0.6 km, geomet¬ 
ric albedo 0.065±0.010 and is a probable C-type (Sheppard 
and Trujillo 2015). The nearly circular orbit at 3.146 AU 
is consistent with membership of the 2 to 3 Gyr old Hy- 
giea family. While the origin of the activity is unknown, 
the 3.33 hr rotation period and 0.45 magnitude lightcurve 


Fig. 8.— Image of 288P (top) and corresponding syndyne plot 
(bottom). A syndyne is the locus of positions of particles of one 
size released from the nucleus at zero speed over a range of times. 
From Hsieh et al. (2012b). 

range suggest, by analogy with 133P (Jewitt et al, 2014a), 
that rotation may play a role. 

2.11. P/2013 R3 (Catalina-PANSTARRS), Tj = 3.185 

P/2013 R3 is a dust enshrouded outer-belt asteroid ob¬ 
served in a state of disintegration (Jewitt et al, 2014b, 
c.f. Figure |^. Ten distinct components were detected in 
the interval from 2014 October to December, with a veloc¬ 
ity dispersion between fragments of order 0.3 to 0.5 m s“^. 
Because of dust contamination in 2014 it is only possible 
to set upper limits to the size of the fragments. Assuming 
geometric albedo 0.05, the four largest have radii < 200 m. 
The gravitational escape velocities of the largest fragments 
(assuming density p = 1000 to 3000 kg m“®) are 0.15 to 
0.25 m s“^, comparable to the measured fragment velocity 
dispersion. The fragmentation took place successively over 
a period of several months at least, excluding an impact as 
the cause. Disruption due to the pressure of a sub-surface 
volatile reservoir was also excluded, although the contin¬ 
ued dust activity of the fragments may have been caused by 
the sublimation of newly exposed ice. Rotational breakup 
seems the most likely cause of the catastrophic disruption 
of R3, supported by the close-to-escape speed relative ve¬ 
locities of the fragments (Hirabayashi et al, 2014b). 

2.12. 133P/(7968) Elst-Pizarro, Tj = 3.184 

133P was first observed to be active in 1996, exhibiting 
a long, narrow dust trail with no visible coma, making it the 
first known (and currently best-characterized) active aster¬ 
oid in the main asteroid belt. At first suspected to be the 
product of a collision (Boehnhardt et al, 1998, Toth 2000), 
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below the ~2 m s“^ gravitational escape speed (e.g. mm 
sized particles have v ^ 6 cm s“^). This is inconsistent 
with a pure sublimation origin, and appears to require cen¬ 
tripetal assistance and a spatially subdivided sublimating 
surface {Jewitt et ah, 2014a). 

The nucleus has a rotation period of 3.471 ± 0.001 hr 
{Hsieh et al, 2004), a spectrum similar to those of B- or 
F-type asteroids (e.g., Bagnulo et al., 2010; Licandro et al, 
2011a), and an optical albedo of pji ^ py ~ 0.05 {Hsieh et 
al., 2009; Bauer et al., 2012). 



Fig. 9.— Evolution of the components of P/2013 R3 over six 
weeks in 2014. The components, labeled A, B and C, show pro¬ 
gressive splitting and are enveloped in a 10* kg debris cloud. Scale 
bars mark 5000 km at the distance of the object. Adapted from Je¬ 
witt et al, (2014b). 


2.13. 176P/(118401) LINEAR (118401), Tj = 3.167 

176P/LINEAR was observed to exhibit a short fan¬ 
shaped tail over a single, month-long interval in 2005 
{Hsieh et al., 2011a). During this time, the object was 
about 30% brighter than the bare nucleus, leading to an im¬ 
plied dust mass ~10® kg. The properties of the dust can 
be approximately matched by models in which the char¬ 
acteristic particle size is 10 pm, the ejection speed ^5 m 
s“^ and the dust production rate ~0.07 kg s“^, all simi¬ 
lar to values inferred in 133P. Activity was sought but not 
detected in 2011 {Hsieh et al, 2014) and thee water produc¬ 
tion rate was spectroscopically limited to Qh^O < 4x 10^^ 
s“^ (< 1.0 kg s“^) in Herschel observations 40 days after 
perihelion on UT 2011 August 8 {de Val-Borro 2012). The 
4.0±0.4 km diameter nucleus rotates with a period near 
22.2 hi- (Table|^. 

2.14. 238P/Read, Tj = 3.152 


it has since been observed to be active on three additional 
occasions in 2002, 2007, and 2013, and inactive inbetween) 
(e.g., Hsieh et al, 2004, 2010; Jewitt et al, 2014a). 



Fig. 10.— 133P imaged UT 2002 September 07. From Hsieh et 
al, (2004). 


The second active asteroid in the main belt, 238P, exhib¬ 
ited a strong coma and dust tail when discovered in October 
2005 at a heliocentric distance of i? ^ 2.4 AU (Figure [TT]i. 
Like 259P, the nucleus of 238P is tiny, with a diameter of 
~0.8 km {Hsieh et al., 201 lb). It was observed to be active 
in both 2005 and 2010, with a period of inactivity in be¬ 
tween, with a coma dust mass of order 10® kg and a produc¬ 
tion rate estimated (from published photometry) near ^0.1 
kg s“^. Also like 259P, 238P is dynamically unstable, with 
a survival time on the order of 20 Myr, although unlike 259P 
which may have been recently implanted at its current loca¬ 
tion, 238P is hypothesized to have diffused in eccentricity 
from its original location within the Themis family (which 
also contains 133P, 176P, and 288P) due to its proximity to 
the 2; 1 mean-motion resonance with Jupiter (Haghighipour 
2009). 


The appearance of 133P when active is typically that of a 
point-like nucleus with a thin tail (or “trail”) of dust follow¬ 
ing in the projected orbit (Figure [T0| ). Order of magnitude 
dust mass loss rates, inferred from surface photometry of 
the tail {Hsieh et al., 2004) reach ~0.02 kg s“^, while com¬ 
parable upper limits in gas are inferred from spectroscopy 
{Licandro etal., 2011a). The thin tail indicates that particles 
are ejected very slowly from the nucleus with characteristic 
speeds ^1.8 m s“^, where is the grain radius 

expressed in microns {Jewitt et al., 2014a). Curiously, par¬ 
ticles larger than a few microns leave the nucleus at speeds 


2.15. P/2012 T1 (PANSTARRS), Tj = 3.134 

P/2012 T1 was discovered to exhibit a diffuse coma and 
a featureless fan-shaped antisolar tail in October 2012. Pho¬ 
tometric monitoring showed that the total scattering cross- 
section of the comet’s coma and tail doubled over a period 
of about a month, remained approximately constant for an¬ 
other 1.5 months, and then declined by ~60% over the next 
1.5 months {Hsieh et al., 2013). Spectroscopic observations 
using the Keck I telescope found an upper limit CN produc¬ 
tion rate of Qcn < 1-5 x 10^® mol s“^, and no evidence of 
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Fig. 11.— 238P imaged on UT 2005 November 10. From Hsieh 
and Jewitt (2006). 


absorption at 0.7 /xm that would indicate the presence of hy¬ 
drated minerals {Hsieh et al, 2013), while Herschel Space 
Telescope observations were used to set an upper limit H 2 O 
production rate of QH 2 O < 7.6 x 10^® mol s ^ {O’Rourke 
et al, 2013). Dust modeling by Moreno et al. (2013) indi¬ 
cated that dust production began near perihelion and lasted 
for a period of ~ 4 — 6 months, with a total ejected dust 
mass on the order of 10^ kg for maximum grain sizes of 
a = 1 — 10 cm. 

2.16. 313P/Gibbs (2014 S4), Tj = 3.132 

Discovered UT 2014 September 24 this object has semi¬ 
major axis 3.156 AU and is located near at least eight other 
active asteroids in the outer belt (Figure |^. The nucleus 
radius is 500 m (albedo 0.04 assumed). Prediscovery ob¬ 
servations from 2003 (Figure [T^ reveal 313P as only the 
third object, after 133P and 238P, to be active in more than 
one orbit, consistent with mass loss driven by the sublima¬ 
tion of ice. It displays a fan-shaped dust tail in both 2003 
and 2014 that is well approximated by syndyne dust emis¬ 
sion models, indicating that the dust is ejected over a period 
of at least 3 months during each active episode {Jewitt et al, 
2015, Hsieh et al, 2015b, Hui and Jewitt, 2015). The ob¬ 
ject is found near two three-body mean-motion resonances 
with Jupiter and Saturn (11J-1S-5A and 10JH-12S-7A), and 
its orbit has been found to be intrinsically chaotic with a 
Lyapunov time of Ti = 12 000 yr, yet numerical simula¬ 
tions show that it is stable over at least 50 Myr. 313P is 
the second active asteroid, after P/2012 Tl, to be associated 
with the ~155 Myr-old Lixiaohua asteroid family {Hsieh et 
al, 2015b). 

2.17. P/2010 R2 (La Sagra), Tj = 3.098 

The object was observed to be active from 2010 Septem¬ 
ber to 2011 January, at i? = 2.6 to 2.7 AU, and has a mea¬ 
sured nucleus size of r = 0.55 ± 0.05 km {Hsieh, 2012c). 
Moreno et al, (201 lb) inferred dust production at the peak 
rate of ~4 kg s“^, with centimeter-sized particles ejected 
at about 0.1 to 0.2 m s“^. A limit to the outgassing rate 
Qcn <3x10^^ s“^ (corresponding to ^3 kg s“^ in wa¬ 
ter) was placed spectroscopically {Hsieh et al, 2012c). The 
latter authors found (neglecting possible non-gravitational 
forces due to outgassing) that the orbit of P/2010 R2 is sta¬ 



Fig. 12.— Prediscovery image of 313P when at 2.47 AU on UT 
2003 October 23. White and grey arrows show the directions of 
the projected anti-solar and heliocentric velocity vectors, respec¬ 
tively. From Hui and Jewitt (2015). 


ble on timescales ^100 Myr and argue that this object was 
likely formed in-situ. 

2.18. 107P/(4015) Wilson-Harrington, T, = 3.083 

This object showed a prominent diffuse tail about 2' in 
length on a blue sensitive photographic plate taken 1949 
November 19, when at i? = 1.148 AU (Fernandez et al, 
(1997). A red-sensitive plate taken nearly simultaneously 
shows only a hint of this tail. 107P appears strongly trailed 
on both plates owing to its non-sidereal motion. The tail 
color {B — R = —1) is too blue to be caused by scattering 
from dust, and the tail position angle (~ 15° from radial 
to the sun) is also inconsistent with the expected direction 
of a dust tail blown by radiation pressure. 107P was re¬ 
observed on November 22 and 25 but then showed no trace 
of a tail {Cunningham, 1950), and no comet-like activity 
has been reported since {Chamberlin et al, 1996, Ishiguro 
et al, 2011c). We regard 107P as the least convincing ex¬ 
ample of the active asteroids, both because the observations 
are old and unrepeated and because 107P lies very near the 
Tj > 3.08 cut-off Bottke et al, (2002) used a statistical 
dynamical model to conclude that there is a 4% chance that 
107P is a captured Jupiter family comet. 

2.19. Related Observations 

Features in the 3.0 to 3.5 pm reflection spectra of large 
main-belt asteroids (24) Themis (diameter 198±20 km) and 
(65) Cybele (273±12 km) have been interpreted as absorp¬ 
tions due to O-H stretch in water ice and C-H stretch in an 
unidentified organic molecule {Campins et al, 2010, Rivkin 
et al, 2010, Licandro et al, 2011b). The shape of the band 
in Themis requires a thin, widespread or global ice film 
(“frost”) only lOOA to lOOOA thick {Rivkin et al, 2010). 
However, such a thin film would be highly unstable to sub¬ 
limation and lead to gas production rates 10® to 10® kg s“^, 
violating observational limits set spectroscopically (<400 
kg s“^, Jewitt and Guilbert-Lepoutre 2012). The latter au¬ 
thors showed that, if water ice exists on these asteroids, it 
must be of high albedo (>0.3) and spatially confined to re- 
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gions far from the subsolar point (e.g. near the poles on 
objects having small obliquity) in order to keep sublima¬ 
tion below observational limits. The spectra of Themis and 
Cybele are similar to that of the iron-rich mineral goethite 
(Beck et al, 2011) but this offers a less plausible expla¬ 
nation because goethite is rare in meteorites and normally 
attributed to weathering of the meteorites once fallen on the 
Earth. 

3. MECHANISMS 

In this section we discuss possible mechanisms able to 
eject dust from asteroids. The group properties offer few 
clues, given the modest size of the active asteroids sample 
(c.f Table Licandro et al, (2011a), for example, sug¬ 
gest that the color distribution of active asteroids is different 
from that of classical cometary nuclei. Even if future mea¬ 
surements prove this to be true, its significance is unclear 
given that a wide range of mechanisms are known to drive 
the observed activity in these bodies (Jewitt 2012). Eurther- 
more, while they are discussed separately for clarity, it is 
likely that different mechanisms operate together in real ob¬ 
jects. (Eor example, the loss of particles produced in (3200) 
Phaethon by thermal disintegration is likely assisted by ro¬ 
tation and radiation pressure sweeping. Sublimation driven 
mass loss from 133P is probably also rotation assisted). 

3.1. Rotational Mass Loss 

Eor a sphere of density p the critical period at which the 
gravitational acceleration equals the centripetal acceleration 
at the equator is 


Pc = 



( 2 ) 


where G is the gravitational constant. Eor example, with 
p = 1000 kg m“^. Pc = 3.3 hr. The critical period is in¬ 
dependent of the asteroid size, is shorter for higher densi¬ 
ties and longer for elongated spheroids in rotation about a 
minor axis (by a factor approximately equal to the ratio of 
the long to short axes of the body). The rotation and the 
shape of a fluid body are related, in equilibrium, by the 
classic MacLaurin and Jacobi ellipsoid series. However, 
while most asteroids may have been severely weakened by 
repeated impact fracturing and the formation of a rubble- 
pile structure, they are not strengthless, van der Waals and 
other weak forces can imbue a rubble pile with a cohesive 
strength while friction can provide resistance to deforma¬ 
tion even in the absence of cohesion, and critical periods 
much shorter than given by Equation are possible (Hol- 
sapple 2007, Sanchez and Scheeres 2014). Calculation of 
the rotational stability of rubble-piles remains an important 
but challenging and under-constrained problem in asteroid 
science. 

Evidence for rotational instability comes from the obser¬ 
vation that most asteroids larger than a few hundred meters 
in size rotate more slowly than a “barrier” period at about 
2.2 hr (Warner et al, 2009). This is widely interpreted as 


meaning that faster rotating asteroids have lost mass or even 
disrupted due to centripetal forces. Recent observations of 
active asteroids 31 IP and P/2013 R3 (and possibly P/2010 
A2) appear to show asteroids losing mass rotationally, so 
that we can begin to study this process as it happens. 

In principle, the rotation rates of asteroids can be driven 
to critical values by external torques exerted by the grav¬ 
ity of other objects, by chance impacts, by outgassing and 
by electromagnetic radiation. In the main belt, gravitational 
torques are negligible (except within binary and other mul¬ 
tiple systems). Collisions add angular momentum stochasti¬ 
cally, leading to a slow random walk towards larger angular 
momenta. Outgassing from icy asteroids can be very effi¬ 
cient in changing the spin but, for ice-free asteroids, radia¬ 
tion torques offer the most potential for driving the rotation 
steadily up to critical values. 

Power absorbed from the Sun by the surface of an aster¬ 
oid is re-radiated to space as heat. Asteroids are aspheri- 
cal and anisothermal, causing the radiated infrared photons 
to be emitted anisotropically. The angle-averaged momen¬ 
tum carried by thermal photons per second corresponds to a 
net reaction thrust on the asteroid, known as the Yarkovsky 
force. It has important dynamical consequences on small 
main-belt asteroids. If the vector representing the net force 
does not pass through the center of mass of the asteroid, the 
result is the so-called YORP torque, which can change the 
magnitude of the spin and excite precession (see the chapter 
by Vokrouhlicky et al in this volume). 

We focus on changes in the magnitude of the spin, and 
estimate the relevant timescale from the ratio of the rota¬ 
tional angular momentum, L, to the torque, T. Ignoring the 
vector nature of these quantities, for simplicity, the torque 
is proportional to the number of photons radiated per sec¬ 
ond, which varies in proportion to where r„ is the 

asteroid radius and R is the distance from the Sun. Torque 
also depends on the moment-arm, defined as the perpendic¬ 
ular distance between the instantaneous direction of the net 
force and the center of mass. Statistically, at least, we ex¬ 
pect the moment arm oc r„. Together, these dependences 
give T oc r^jR^. Meanwhile, the spin angular momentum 
is L oc Mr'^uj, where M is the body mass and uj the an¬ 
gular rotation rate, related to the rotational period, P, by 
w = 27r/P. Substituting M oc pr^ gives L oc pr^co. Ei- 
nally, the timescale for YORP to change the angular mo¬ 
mentum is Ty = L/T, or Ty = Kpr'^R^ui, where AT is a 
constant. 

The value of constant K depends sensitively on the body 
shape, surface texture, thermal properties and spin vector of 
the asteroid. We estimate its magnitude from measurements 
of YORP acceleration in five asteroids (Rozitis and Green 
2013, Lowry et al, 2014), scaled to assumed density p = 
2000 kg m“^ and spin period P = 5 hr, to obtain 


Ty{Myr) ~ 4 


■ rn ■ 

2 

R 

. 1 km- 


3AU 


(3) 


Order of magnitude deviations from Equation Q can re- 
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suit from peculiarities of the asteroid shape, texture, ther¬ 
mal properties and spin. Still, Equation Q gives a crude 
but useful estimate, by showing that the angular momentum 
of a 1 km radius asteroid at 3 AU can be modified signifi¬ 
cantly on timescales of just a few Myr. A sustained YORP 
torque can drive a kilometer-scale body to rotational insta¬ 
bility in a few times Ty, and the spins of sub-kilometer, low- 
density, near-Sun asteroids should be particularly suscepti¬ 
ble to YORP torques. Setting Ty < 4.5 x 10® yr in Equation 
([^ shows that asteroids as large as r„ ^ 30km can be af¬ 
fected by YORP in the age of the solar system, provided 
they survive against collisions for such a long time. Most 
of the objects in Table Q are small enough to be potentially 
affected by YORP. 

A major complication is that YORP torques can change 
not just the rotation state but the shape of an asteroid, as ma¬ 
terial slides, bounces or rolls towards the rotational equator 
on its way to escape (e.g. Harris et al, 2009, Statler 2009). 
Changing the shape affects the magnitude and possibly the 
direction of the YORP torque (Cotto-Figueroa et al, 2014), 
creating a feedback loop that can change the timescale for 
spin-up relative to the value in Equation ([^. Evidence for 
equatorial accumulation of particulate matter is present in 
images of UEO-shaped Saturnian satellites (c.f. Eigure [T3|) 
and in the shapes of some rapidly-rotating near-Earth aster¬ 
oids, as determined by radar {Ostro et al, 2006). 

In part because of the complex interplay between torque, 
spin, material properties, shape and mass loss, no compre¬ 
hensive models of asteroid evolution and disruption by ro¬ 
tation have been computed, although many interesting ap¬ 
proximations exist {Holsapple 2007, 2010, Walsh et al, 
2008, 2012, Jacobson and Scheeres 2011, Marzari et al, 

2011, Hirabayashi et al, 2014a, Sanchez and Scheeres 

2012, 2014, Cotto-Figueroa et al, 2014, Scheeres 2015). 
These models, not always in agreement with one another, 
attempt to investigate the role of small torques in the pro¬ 
duction of distinctive body spins and shapes, in mass- 
shedding and structural failure (c.f. Eigure [T^, and in the 
formation of both binaries and unbound pairs (Pravec et al, 
2010, Moskovitz 2012, Polishook 2014). 

The cohesive strength is a key parameter in all such mod¬ 
els because, given sufficient cohesive strength, an asteroid 
can resist rotational forces at an arbitrarily high rotation rate 
(e.g. Rozitis et al, 2014). Eortunately, new observations of 
active asteroids allow us to obtain a useful estimate of the 
cohesive strength. 

Eirst, we consider a toy model of two bodies in contact, 
both of density p and having dimensions rp and r^, with 
Cp ^ Tg. Assume that the smaller body sits on the rota¬ 
tional equator of the larger body and take the area of con¬ 
tact between them to be ~ r^. Let the critical frequency 
of rotation at which gravitational attraction is balanced by 
centripetal acceleration be ujc- Then, when rotating with an 
angular frequency w > ojc, the minimum cohesive strength 
needed to bind the two bodies together may be written 



Eig. 13.— Two views of UFO-shaped Saturnian satellite Atlas 
showing a dust skirt around the equator where the effective grav¬ 
itational attraction to the body approaches zero. Atlas is a slow- 
rotator (period 14 hr), hut it fills its Roche lohe within the Saturn- 
inan system, making it a suitable analog for a body at the stability 
limit {Minton 2008). The paucity of craters indicates youth, which 
results from the mobility of surface dust. A possible landslide is 
visible on the dust skirt; corresponding landslides or avalanches 
leading to escape are likely responsible for the impulsive tail ejec¬ 
tions in 31 IP (c.f. Figure |^). The bulk density ~400 kg m“^. 
Credit: NASA/JPL/SSI. 


S ~ ~ (4) 

^ s 

where rris is the mass of the smaller component. If de- 
tachment occurs when >> the smaller mass uig 
will leave with a relative velocity, AU, comparable to the 
instantaneous equatorial velocity of the primary, namely 
AU = TpUj. Substituting for AU and rUg = pr^ into Equa¬ 
tion 0 we obtain 

^-p(^^^(AU)2 (5) 

The six measured components of P/2013 R3 (Table 2 of 
Jewitt et al, 2014b) are equivalent in volume to a single 
precursor body of dimension rp = 350 m. Eour of these six 
nuclei have Vg ^ 200 m and are separating at characteris¬ 
tic speeds 0.2 < AU < 0.5 m s“^. Assuming density p 
= 1000 kg m“^ Equation 0 gives an estimate of the co¬ 
hesive strength as 16 < S' < 140 N m“®. A mathemati¬ 
cally more refined analysis, in which the initial body shape 
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Structural failure Surface shedding 



Fig. 14.— Schematic diagram distinguishing (left panel) struc¬ 
tural failure, as may be seen in P/2013 R3, from (right panel) sur¬ 
face shedding (for which 3IIP is a candidate). From Hirabayashi 
and Scheeres, (2014). 

is represented by a rotational ellipsoid, gives 40 < S ^210 
N m~^ (Hirabayashi et al, 2014b). In both cases, these 
very small cohesive strengths are comparable to the ~25 N 
m~^ strength expected of a rubble pile bound by van der 
Waals forces acting on the smallest particles (Sanchez and 
Scheeres 2014). The small size of the P/2013 R3 precursor 
body is compatible with a short YORP timescale (Ty ~ 0.5 
Myr by Equation]^, lending credibility to the YORP spin- 
up hypothesis for this object. Active asteroid 31 IP has also 
been interpreted (Jewitt et al, 2013c,2014c) as rotational 
instability of a different sort, in which only surface particu¬ 
lates escape from the rotating central body (Hirabayashi et 
al, 2014). 



Diameter [km] 


a major disruptive role. Indeed, the rotational disruption 
rate may exceed the impact destruction rate for asteroids 
<1 km in diameter (Figure [T5]l. However, many questions 
remain about the process. In 31 IP, what determines the in¬ 
tervals between successive dust releases and for how long 
will they continue? More specifically, does the distribution 
of dust releases resemble the power law spectra of self- 
organized critical sandpiles (Laurson et al, 2005), as might 
be expected if avalanches are responsible? Could this ob¬ 
ject be in the process of forming a binary (c.f. Walsh et 
al, 2008, Jacobson and Scheeres 2011) or has it already 
formed one? Will the nucleus eventually be driven to struc¬ 
tural instability, as has occurred in P/2013 R3? Or could 
31 IP be a fragment of a rotationally disrupted precursor 
body whose components are so widely spread as to have 
escaped notice? 

3.2. Impacts 

The average velocity dispersion among main-belt aster¬ 
oids is t/ ~ 5 km s“^, enough to cause vaporization as well 
as shattering and, potentially, disruption upon impact. Ex¬ 
periments at laboratory scales show that the ratio of nie, the 
mass ejected faster than speed v, to the projectile mass, M, 
is 



with A ~ 0.01 and, very roughly, a ~ -1.5 (Housen and 
Holsapple 2011). For an assumed spherical projectile, the 
mass is M = Anpr^/i, where rp is the radius. Equation (j^ 
is valid provided the impact does not have enough energy 
to completely disrupt the target asteroid. For example, an 
impact at 5 km s“^ into a 1 km radius body, with escape 
velocity Ue ~ 1 m s“^, would have an ejecta to projectile 
mass ratio rrie/M ^ 3500 by Equation (j^. The mass of 
ejecta and the scattering cross-section, C, are related by 

We = XPaC (1) 

where d is the weighted mean particle size in the ejecta and 
X is a dimensionless constant of order unity. Very modest 
projectiles are capable of creating ejecta with a substantial 
cross-section. We write C = where is the radius 

of the target asteroid and / is a constant. The case / = 1 cor¬ 
responds to ejecta with a cross-section equal to that of the 
target asteroid. Such an impact would result in a doubling 
of the asteroid brightness, and therefore should be readily 
detected. We set v = (STrCp/Sj^/^rn, the gravitational es¬ 
cape speed from the target asteroid, and combine Equations 
(|^ and (0 to obtain 


Fig. 15.— Model rates of asteroid disruption due to impact 
and rotational bursting as a function of asteroid diameter. From 
Marzari et al, (2011). 


Tp = 




( 8 ) 


This first determination of the cohesive strength, if 
broadly applicable, shows that rotation of asteroids can play 
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Equation (|^ assumes that the densities of the target and pro¬ 
jectile are equal. Substituting the above parameters. Equa¬ 
tion (|^ reduces to 




Tp = 0.29 (9) 

Equation is plotted in Figure ( |16[ ) for / = 0.1, 1, 10 
and 100, corresponding to changes in the apparent magni¬ 
tude due to impact of Am ~ 2.51og2o/ = 0.1, 0.7, 2.5 
and 5 magnitudes, respectively. The Figure shows that a 
20 cm projectile striking a 1 kilometer asteroid would eject 
enough material to double the total cross-section and bright¬ 
ness. Very modest impacts should create observable ejecta 
signatures. 



Fig. 16.— Projectile radius (y-axis) needed to generate ejecta 
with cross-section / = 0.1, 1, 10 and 100 times the geometric 
cross-section of the target asteroid, whose radius is plotted on the 
x-axis (c.f. Equation 1^. Impacted asteroid (596) Scheila (vp = 
65 km) brightened by about 1 magnitude following impact. The 
figure shows that the projectile radius was several lO’s of meters. 

The number of boulders in the asteroid belt is very uncer¬ 
tain, because these small objects are too faint to be observed 
directly from Earth. Models extrapolating from larger sizes 
suggest that the number of objects with rp > 1 m is roughly 
10i3±i (^Q’Brien and Greenberg, 2005). The rates at which 
small-scale impacts occur is correspondingly uncertain by 
at least an order of magnitude. 

Survey observations by Denneau et ai, (2014) have been 
interpreted, through a simple model, to show a discrepancy 
with the rate of asteroid disruption predicted by Bottke et 
ai, (2005). 


3.3. Thermal Disintegration 

Forces caused by thermal expansion can exceed the frac¬ 
ture strength of a material, leading to cracking and the pro¬ 
duction of dust. On a low gravity body, excess thermal 
strain energy can eject fragments faster than the gravita¬ 
tional escape speed, leading to a net loss of material {Je- 
witt 2012). In a uniform material, stresses arise from dif¬ 
ferential thermal expansion across a temperature gradient. 
Most minerals are non-uniform in structure, consisting of 
compositionally distinct grains. Non-uniform materials are 
susceptible to fracture because of differential expansion be¬ 
tween the component materials, with or without a tempera¬ 
ture gradient. The thermal skin depth in the material (across 
which temperature gradients are particularly marked, and 
variable) is related to the period of the temperature varia¬ 
tion, T, and to the thermal diffusivity, k, by f ~ 

For example, an asteroid with diffusivity k = 10“® m^ s“^, 
rotating with period r = 5 hr would have a skin depth £ ^ 
0.1 m. Thermal fracture occurring within a thermal skin 
depth of the physical surface may contribute to the produc¬ 
tion of regolith (Delbo et ai, 2014). Desiccation (the loss of 
bound water from a hydrated mineral) is a separate process 
leading to fracture and dust production when the stresses in¬ 
duced by shrinkage exceed the material strength. On Earth, 
mudcracks in dry lake beds made of clay minerals are the 
most easily recognized examples of desiccation cracking. 

Among the active asteroids, the most convincing evi¬ 
dence for thermal disintegration is from (3200) Phaethon 
{Li and Jewitt 2013, Jewitt and Li 2010, 2013). The ec¬ 
centric orbit and high obliquity (Phaethon’s pole is at eclip¬ 
tic coordinates A,/3 = 85°±14°, -20°±10°, Ansdell et al., 
2014) lead to highly non-uniform surface heating, with peak 
temperatures near 1000 K at perihelion {Ohtsuka et ai, 
2009). The temperature is further modulated by Phaethon’s 
rapid rotation (3.6 hr period), with a diurnal variation of 
hundreds of Kelvin. Furthermore, the escape speed from 
the ^5 km diameter nucleus is a few meters per second, 
sufficiently small that grains produced by fracture can eas¬ 
ily be ejected {Jewitt 2012). While Phaethon’s composition 
is unknown, hydrated mineralogies (montmorillonite clays) 
have been suggested based on the blue optical/near-infrared 
reflection spectrum {Licandro et ai, 2007). Hydrated mate¬ 
rials, if present, would be highly susceptible to desiccation, 
shrinkage and fracture. Nine asteroids with perihelia <0.25 
AU (subsolar temperatures >800 K) have been searched for 
activity, setting upper limits to the mass loss rate from <0.1 
to 1 kg s“^ {Jewitt 2013). 

On a larger scale, thermal disintegration may be im¬ 
portant in dust disks surrounding white dwarf stars. Such 
disks are inferred from photospheric, heavy-element pollu¬ 
tion {Jura and Xu 2013) 

3.4. Sublimation of Ice 

Ice exposed at the surface of an asteroid at distance R 
AU from the Sun will sublimate, in equilibrium, at a specific 
rate dm/dt (kg m~^ s“^), given by the solution to 
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g cc»(t)) = ear *UT) — 


( 10 ) 


where the term on the left represents the absorbed solar 
power per unit area, the hrst term on the right represents 
thermal radiation to space per unit area and the second term 
represents the power used per unit area in breaking bonds 
to sublimate ice. Fq = 1360 W m ^ is the solar constant, 
A is the Bond albedo, e is the emissivity of the asteroid 
and L{T) is the latent heat of sublimation of ice at temper¬ 
ature T. Here 9 is the angle between the Sun and the sur¬ 
face normal as seen from the sublimating surface. Values of 
cos{9) range from 1/4 (the isothermal case, corresponding 
to the lowest possible temperatures) to 1 (for normal illu¬ 
mination at the sub-solar point, the maximum temperature 
case). Conduction is ignored for simplicity - its effect will 
be to reduce dm/dt relative to the value computed from 
Equation ( [T0| i by an amount that depends on the conductiv¬ 
ity (measurements of small body regoliths suggest that the 
conductivity is very small, so that the error incurred by its 
neglect should be minor). Equation ( [T0) i can be solved in 
combination with the Clausius-Clapeyron relation for wa¬ 
ter ice, to evaluate dm/dt{R), as in Eigure (17 1 . The Eig- 
ure also shows (on the right hand axis) the rate of reces¬ 
sion of the sublimating surface due to mass loss, given by 
dr/dt = p~^dm/dt (m s“^). 

Maximum sublimation rates vary by nearly an order of 
magnitude from ^10“^ kg m“^ s“^ at 2 AU to ^10“® kg 
m“^ s“^ at 3.5 AU, corresponding to the inner and outer 
edges of the main -belt. The corresponding surface re¬ 
cession rates are from a few meters to a few x 0.1 meters 
per year. Water ice is thermodynamically unstable when 
exposed to sunlight at asteroid belt distances, but can be 
preserved on billion-year timescales beneath meter-thick or 
greater layers of porous refractory debris (Fanale and Sal- 
vail 1989, Schorghofer 2008). Sublimation at rates suffi¬ 
cient to launch dust from a small body would then require 
that the protective regolith be removed, perhaps by a small 
impact or some other disturbance of the surface. In this sce¬ 
nario, sublimation would proceed from an exposed region 
until a new refractory covering, consisting of particles too 
large to be lifted by gas drag (a so-called “rubble mantle”), 
chokes off the flow. The ratio of the lifetime of the sublimat¬ 
ing surface patch to the interval between successive impacts 
(or other surface disturbances) defines the “duty cycle”, 

Eor example, the timescale for formation of a rubble mantle 
on a kilometer-scale body at 3 AU is estimated from a sim¬ 
ple model (Eigure 4 of Jewitt 2002) at ^1 yr, increasing to 
10 yr at 4 AU. The interval between impacts of meter-scale 
projectiles (needed to expose ice sufficient to account for 
the mass loss) onto this body is perhaps t ^ 10^ yr, giving 
fd ~ 10“^ to 10“^. Neither the re-sealing timescale nor 
the excavation timescale can be specihed with conhdence, 
but /d <C 1 is assured. As noted in Section [T2l impacts at 5 
km s“^ are highly erosive, with an ejecta to projectile mass 
ratio rrie/M ~ 3500 for a 1 km target asteroid (Equation 
1^. At this rate, the asteroid can survive for ~10® yr. 



^ 10 ' 


^ 10 ” 


^ 10 ” 


10 ” 


10 ” 


10 ” 


c 

a 




Heliocentric Distance, R [AU] 


Eig. 17.— Solutions to Equation 1 10' for the extreme cases of 
subsolar sublimation (6^ = 0°) and for isothermal sublimation. The 
right hand axis shows the rate of recession of the sublimating sur¬ 
face due to mass loss, given by dr/dt = p~^ dm/dt (m s“^), 
where p is the density, here taken to be 1000 kg m~®. From Jewitt 
( 2012 ). 


3.5. Radiation Pressure Sweeping 

Particles lifted above the surface of an asteroid are sus¬ 
ceptible to solar radiation pressure and, if sufficiently small, 
can be blown away. The critical size for radiation pressure 
sweeping on a non-rotating asteroid is {Jewitt 2012) 


a0{pLm) = 10 


1 km 



( 11 ) 


At i? = 3 AU, a 1 pm radius grain detached from the surface 
(by any process) can be blown away from a non-rotating 
asteroid 1 km in radius. 


3.6. Electrostatics and Gardening 

Dielectric surfaces exposed to ionizing (UV and X-ray) 
solar radiation develop a positive electric potential (of per¬ 
haps 5 to 10 V) through the loss of photoelectrons. Con¬ 
versely, electrons impacting in shadowed regions imbue a 
locally negative potential. The resulting electric fields near 
shadow boundaries can grow large enough (~10 V m~^ to 
100 V m“^) to mobilize dust. On the Moon, horizon glow 
(Rennilson and Criswell 1974) and dust impact counter ex¬ 
periments {Berg et al, 1976) show that charging effects 
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near the terminator lift 10 /rm sized dust particles to me¬ 
ter heights and greater, at implied ejection speeds, ~1 m 
s“^. Smooth dust ponds on asteroids strongly suggest elec¬ 
trostatic mobilization (but not ejection) of grains in a much 
lower gravity environment (Poppe et ai, 2012, see also the 
chapter by Murdoch et al, in this volume). 

On small asteroids, electrostatic forces can potentially 
accelerate particles to speeds exceeding the gravitational es¬ 
cape speed. To within a numerical factor of order unity, 
particles smaller than a critical radius 

ae{pm) ~ 2 (12) 

can be electrostatically ejected against gravity from a non¬ 
rotating body of radius {Jewitt 2012). 

Larger particles are too heavy to lift against gravity, at 
least on a non-rotating body, whereas smaller particles are 
trapped by cohesive forces on the surface {Hartzell and 
Scheeres 2011). Electrostatically ejected particles should 
therefore be confined to a narrow range of sizes near Ug. 
Note that Equation 0 is independent of the heliocentric 
distance because the potential, V, depends on photon en¬ 
ergy, not flux. Since, in the asteroid belt at 2 to 3 AU, 
Oe < 0 / 3 , particles launched electrostatically should be 
picked up by solar radiation pressure and propelled into 
comet-like dust tails. 

Although electrostatic ejection must operate at some 
level, in practice the rate of loss of small particles from 
an asteroid will be limited by the rate of their production, 
presumably by micrometeorrte bombardment. This distin¬ 
guishes the role of electrostatics on the asteroids from the 
circumstance on the Moon, where ejected particles fall back 
to the surface under the high lunar gravity, to be repeatedly 
relaunched at each terminator crossing (Berg et ai, 1976). 

We crudely limit the ejected mass flux from asteroids 
as follows. On the Moon, impact overturn of the regolith 
(“gardening”) reaches to depths Sr <1 m on timescale 
t = 1 Gyr (Gault et ai, 1974). Impact-produced frag¬ 
ments that are largely retained on the Moon by gravity will, 
on the small asteroids, be immediately lost. Therefore, a 
strong upper limit to the electrostatic mass loss from aster¬ 
oids can be obtained by assuming that the gardening rates 
are comparable and that none of the gardened layer is re¬ 
tained. The mass loss rate from a spherical asteroid is 
dM/dt ~ Airr'^pSr/t, where is the radius and p is the 
density. Eor an = 1 km asteroid with p = 2000 kg m“^ 
and 5r = 1 m, this gives dM/dt < 10“® kg s“^. This 
is 4 to 6 orders of magnitude smaller than the mass loss 
rates inferred from observations of active asteroids (Table 
1^. Of course, we cannot exclude spatially or temporally 
enhanced micrometeorite impact fluxes in the asteroid belt, 
but enhancements sufficient to bring the mass loss rates up 
to those in Table Q seem unlikely. The only observation 
possibly supporting the existence of widespread, very low 
level mass loss is by Sonnett et al, (2011). If this obser¬ 
vation is confirmed, it might be evidence for very low level 
ejection of dust owing to gardening losses, possibly assisted 


by electrostatic effects. 

3.7. Observational Diagnostics 

Diagnosing the cause of mass loss from any particular 
object is difficult, both because of limitations of the data 
and because of uncertainties in models of the mass loss 
mechanisms. It is often easier to rule out possible expla¬ 
nations than to rule them in. Spectroscopic detections of 
water vapor or another sublimation product would provide 
the most solid evidence of sublimation-driven activity, but 
such evidence has only been obtained for Ceres (Kiippers et 
al, 2014). Unsuccessful attempts have been made to detect 
water vapor or CN emission from several other active aster¬ 
oids (e.g., Jewitt et al, 2009; Licandro et al, 2011a, 2013; 
Hsieh et al, 2012a,b,c, 2013; de Val-Borro et al, 2012; 
O’Rourke et al, 2013). Typical derived upper limits to gas 
loss rates are <1 kg s“^. Unfortunately, the significance of 
these non-detections is inconclusive, for two reasons. Eirst, 
CN is a prominent trace species in the spectra of Kuiper belt 
and Oort cloud comets (with H 2 O/CN ~360); it might be 
strongly depleted in icy asteroids. Second, spectra are typ¬ 
ically obtained long after continuum observations have re¬ 
vealed the mass loss - the gas could simply have escaped on 
a timescale short compared to the radiation pressure sweep¬ 
ing time for dust. 

Another potential indicator is the recoil (or “rocket”) ac¬ 
celeration produced on the nucleus by anisotropic sublima¬ 
tion. Specifically, a non-gravitational acceleration, a, im¬ 
plies a mass loss rate 

dM _ aMn 

dt ~ fevo ^ ’ 

where Mn is the mass of the nucleus and vq is the veloc¬ 
ity of the material ejected. Dimensionless parameter fg ac¬ 
counts for the directional pattern of the mass loss, with fg = 
0 and 1 cotTesponding to isotropic and perfectly collimated 
ejection, respectively. In fact, non-gravitational accelera¬ 
tion of 133P at a = 2.5 x 10“^° m s“^, has been reported 
(Chesley et al, 2010), albeit with only 3 (t statistical signif¬ 
icance. To obtain the smallest possible estimate of dM/dt 
from Equation ( [T3] ), we maximize fg = l and use the sound 
speed in H 2 O gas at the blackbody temperature appropri¬ 
ate to 3 AU, namely vq = 500 m s“^. Approximated as a 
sphere of radius 2.2 km and density 1000 kg m“^, the nu¬ 
cleus mass is = 4x10^^ kg, giving dM/dt ^ 20 kg 
s“^. This rate exceeds the dust production rates measured 
in 133P by two to three orders of magnitude, forcing us to 
conclude either that a has another cause or that the reported 
value should be considered only as an upper limit. Nugent 
et al (2012) measured small semi-major axis drift rates in 
42 near-Earth asteroids consistent with the Yarkovsky (ra¬ 
diation) acceleration. An additional 12 asteroids showed 
drift rates too large to be explained by the Yarkovsky effect. 
If confirmed, the accelerations of these asteroids, none of 
which is a known active asteroid, might suggest anisotropic 
mass loss. 
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Less direct inferences can be made in the absence of 
spectroscopic detection of gas. Activity which repeats on 
the orbital timescale, as in 133P and 238P, is naturally ex¬ 
plained only by sublimation. A prolonged period of dust 
emission is also suggestive of sublimation and difficult to 
explain as the result of a simple impact. However, pro¬ 
longed emission may not be unique to sublimation. Evi¬ 
dence from 31 IP shows that mass shedding can continue for 
months, raising the potential for confusion with protracted 
emission caused by sublimation. The size vs. velocity re¬ 
lation for ejected dust particles also differs between mech¬ 
anisms, with gas drag giving v cx where a is the 

dust radius, and rotational instabilities and impact giving a 
flatter dependence (Figure [TSll. 


Particle Radius [m] 



Fig. 18.— Dust ejection velocities as function of the radiation 
pressure coefficient /3 or particle size. The plot is a compilation 
of published models for individual objects. Such models typi¬ 
cally simulate the motion of dust grains ejected from the nucleus 
over a range of times, particle sizes, and velocities. The best¬ 
fitting parameters are determined from comparison with the ob¬ 
served shapes of the dust tails. Solid lines indicate objects with 
protracted, likely sublimation-driven dust emission, dashed lines 
and the shaded area for P/2010 A2 refer to objects with impulsive 
emission. Arrows indicate upper limits, the shaded area for P/2010 
A2 illustrates that particles of all sizes were ejected with the same 
range of velocities. 

Sublimation will preferentially lift small, micron-sized 
particles that are quickly dispersed by radiation pressure. 
Objects typically have a bright coma and tail for a few 
months, as long as the sublimation is on-going, followed 
by a rapid fading and return to the appearance of a point 
source. During rotational breakup, the ejecta sizes are not 
constrained at all, due to the compensation of gravity by 
the centrifugal force. Also impacts are expected to exca¬ 
vate much larger pieces of debris, depending on the escape 
speed from the parent body. Characteristically, small ob¬ 
jects suspected to have undergone an impact or rotational 
breakup display long-lived trails of cm- to m-sized debris 
along their orbits for many years after the event, because 
such large grains are not efficiently removed by radiation 
pressure. (596) Scheila, though clearly impacted, did not 


develop a debris trail because the high escape speed pre¬ 
vented the ejection of large debris. The absence of de¬ 
bris trails in sublimation-driven active asteroids is different 
from comets, likely because cometary sublimation is much 
stronger and able to lift also large debris. 

Numerical modeling is necessary because the appear¬ 
ance of long-lived activity can be produced either by the lin¬ 
gering of large, slow-dissipating dust grains, or the ongoing 
replenishment of small, fast-dissipating dust grains. Un¬ 
der the right conditions, however, dust modeling can break 
this degeneracy and give us insights into the origin of active 
asteroid activity (e.g., Hsieh et al., 2004; Ishiguro et al, 
201 lb; Agarwal et al., 2013). 

Pinpointing cases of rotational instability can be diffi¬ 
cult because there are few clear predictions about the likely 
appearances of such bodies. The repeated dust ejections 
from 31 IP, for example, can be reconciled with a rotational 
shedding instability, but it is quite likely that a range of mor¬ 
phologies and time-dependent emission profiles can be pro¬ 
duced by this mechanism. Significantly, other mechanisms 
seem incompatible with the observations. The multiple ob¬ 
ject P/2013 R3 is self-evidently a result of break-up, with 
rotation as the likely cause. 

Table [3 summarizes our conclusions about the mecha¬ 
nisms driving activity in each object. For many objects, we 
can eliminate several processes but cannot isolate a unique 
cause. This partly reflects the inadequacies of the data but 
is also a true consequence of the nature of the activity. For 
example, in those cases where sublimation is believed to be 
the primary activity driver, an impact (or other disturbance) 
may be needed to excavate buried ice to trigger sublimation, 
and mass loss may be assisted by rapid rotation. We empha¬ 
size that, in many cases, the entries in Table are expected 
to change with the acquisition of new data. 

3.8. Distribution 

The orbital distribution of the active asteroids may pro¬ 
vide a clue as to the activity mechanisms. There is apparent 
clustering of active asteroids in Figure in the outer belt 
near a ^ 3.1 ±0.1 AU. Of the 13 active asteroids orbiting 
in the main-belt, only two (P/2010 A2 and 3IIP) are lo¬ 
cated inside the 3:1 mean-motion resonance with Jupiter (a 
= 2.5 AU), corresponding to ~15% of the total while, for 
the asteroids as a whole, this ratio is ~35%. However, the 
difference between these distributions is not significant. A 
Kolmogorov-Smirnov test applied to compare the semima- 
jor axes of the active asteroids and the first 500 numbered 
asteroids shows that the distributions have a 17% probabil¬ 
ity of being drawn from the same parent distribution and, 
therefore, are consistent. However, when the comparison is 
made between the distribution of the semimajor axes of ac¬ 
tive asteroids and of main-belt asteroids of comparable size, 
the difference instead becomes highly significant. Specifi¬ 
cally, there is <0.1% likelihood that active asteroids and 
main-belt asteroids with absolute magnitude iT > 18 are 
drawn from the same parent population, according to the K- 
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S test. This different result occurs because of observational 
bias in magnitude-limited sky survey data, which strongly 
inhibits the detection of small asteroids in the outer belt. 
While the active asteroids are subject to this same bias, their 
detectability is also influenced by their activity, as is evi¬ 
dent from the fact that many objects (e.g. P/2010 A2, 31 IP, 
P/2013 R3) were unknown before being discovered in an 
active state. We conclude that the relative paucity of active 
asteroids in the inner-belt and the concentration near 3.1 
AU are real, but larger samples from quantitative sky sur¬ 
veys will be important to better assess the biases. Clustering 
could suggest an origin as collisionally produced fragments 
of an ice-rich precursor (Novacovik et ai, 2012, 2014). 

4. DYNAMICS 

Numerical simulations have been conducted to assess 
the long-term dynamical stability of many of the known 
active asteroids (e.g., Jewitt et ai, 2009; Haghighipour, 
2009; Hsieh et al., 2012a,b,c, 2013; Stevenson et ai, 2012). 
This is of particular interest in the context of understanding 
whether ice-bearing main-belt comets are native to the as¬ 
teroid belt or are implanted interlopers from the outer solar 
system. In most cases, while the median dynamical lifetime 
of short-period comets before ejection from the solar sys¬ 
tem or collision with the Sun was found to be ~ 5 x 10® yr 
(Levison & Duncan, 1994), active asteroids in the asteroid 
belt have been found to be stable over timescales of 10® yr 
or longer. 238P and 259P are notable exceptions, however, 
having been found to be unstable on timescales on the order 
of ^ 10^ yr {Jewitt et al, 2009; Haghighipour, 2009). 

Fernandez et al. (2002) forward-integrated the orbits of 
test particles representing the 202 JFCs known at the time 
and several dynamical clones of each comet (using purely 
gravitational computations). None ended up on orbits with 
both low eccentricity and low inclination like that of 133P. 
Comet 503D/Pigott was seen to evolve onto a main-belt or¬ 
bit, but with high inclination. Interestingly, we know now of 
active asteroids in the main belt with high inclinations (e.g., 
259P and P/2010 R2), and so the possibility that these could 
be JFC interlopers cannot be excluded. In the same work, it 
was shown that detaching objects from the JFC population 
was possible with the inclusion of non-gravitational (NG) 
forces in the simulations, but such large forces were needed 
as to be unrealistic. The size dependence of NG forces plays 
a role in this, since NG forces should be far less effective 
for larger objects, while smaller objects may not be able to 
sustain such large NG forces for very long without disinte¬ 
grating or exhausting their volatile content. 

The aforementioned simulations assume the modern-day 
architecture of the solar system. Some dynamical mod¬ 
els of the early solar system suggest that planet migration 
could have resulted in the emplacement of icy outer solar 
system objects in what is now the main asteroid belt (e.g., 
Levison et ai, 2009; Walsh et al, 2011), although proba¬ 
bly not in low inclination orbits like those of many active 
asteroids. Therefore, particularly when assessing the astro- 


biological significance of icy active asteroids as indicators 
of the composition of the inner protosolar disk, the possi¬ 
bility that these objects could have been implanted at early 
times must continue to be considered (see the chapter by 
Morbidelli et al in this volume). 

Several active asteroids are associated with collisional 
asteroid families and clusters. For example, four active as¬ 
teroids (133P, 176P, 238P, and 288P) have been linked to the 
Themis family either currently or in the past {Toth, 2000; 
Hsieh, 2009; Haghighipour, 2009; Hsieh et al, 2012c), 
while P/2010 A2 and 31 IP are associated with the Flora 
family {Moreno et al, 2010; Jewitt et al, 2013c) and 
P/2012 T1 and 313P with the Lixiaohua family {Hsieh et 
al, 2013, 2015b). The dynamical associations are consis¬ 
tent with photometry showing, for example, that the colors 
of 133P, 176P and 238P are close to those of the Themis 
family {Licandro et al, 2011a) while that of 31 IP is con¬ 
sistent with membership in the Flora family {Jewitt et al. 
2013c). However, dynamical associations are unsurpris¬ 
ing, given that many asteroids are family members, and a 
causal relationship between family membership and activ¬ 
ity cannot be assumed. The color relationships could also 
be coincidental, since S-types are common in the inner-belt 
where P/2010 A2 and 3IIP are found, and C-types in the 
outer belt where 133P, 176P, 238P, and 288P orbit. How¬ 
ever, several active objects have also been linked to much 
smaller, young asteroid clusters, including 133P (part of the 
<10 Myr-old Beagle family; Nesvorny et al, 2008), 288P 
(part of a newly-discovered 7.5±0.3 Myr-old cluster; No- 
vakovic et al, 2012), and P/2012 F5 (part of a 1.5±0.1 Myr- 
old cluster; Novakovic et al, 2014). In the case of active 
asteroids exhibiting sublimation-driven activity, member¬ 
ship in a young cluster suggests a natural mechanism by 
which ice could have been preserved over long timescales 
(i.e., within the interior of a larger parent body) and only 
recently exposed to more direct solar heating. In the case 
of objects exhibiting impact-induced activity, an associated 
young cluster could represent both a source of potential im- 
pactors (given that cluster members will share very similar 
orbits) and possibly a consequence of residing in a region of 
the asteroid belt characterized by high collision rates (e.g., 
Novakovic et al, 2014). 

5. SURVEYS AND RATES 

About 10 of the known active asteroids (Table and 
~10® main-belt asteroids {Jedicke and Metcalfe 1998) are 
larger than 1 km. The ratio of these numbers gives an active 
fraction, for all types of activity-driving mechanism, / ^ 
10“®. However, most asteroids fall near the limiting mag¬ 
nitudes of the surveys in which they are discovered and so 
are never effectively searched for dust emission. Therefore, 
/ ~ 10“® sets a strong lower limit to the active fraction, 
since many active objects must go undetected. Recent, ded¬ 
icated surveys have attempted to remedy this situation by 
searching for comae in well-defined sky surveys. Unfor¬ 
tunately, most of these surveys are so small in scale that 
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they detect no active objects (of the known active aster¬ 
oids, only 118401 LINEAR (1999 RE70) was discovered 
as the result of a targeted survey), allowing 3a upper limits, 
/scr, to be determined. These include < 54 per million 
(Gilbert and Wiegert 2009, 2010), fsa- < 6000 per million 
(Sonnett et al. 2011), and /acr < 50 per million (Waszczak 
et al. 2014). These surveys were based on the attempted 
detection of spatially resolved coma. Cikota et al., (2014) 
instead used photometry to search for the brightening pro¬ 
duced by the ejection of dust, finding five candidates from 
a search of ~75 million observations of ~300,000 aster¬ 
oids. These candidates remain unconfirmed and it is not 
clear which method, resolved imaging vs. integrated pho¬ 
tometry, offers the better sensitivity to mass loss. 

Hsieh et al. (2015a) concluded from a sample of 30,000 
objects observed near perihelion (giving two detections) 
that / ~ 100 per million (in the outer belt), which we 
take as the best current estimate of the fraction of aster¬ 
oids which are measurably active at any instant. This is 
still a lower limit to the true active fraction, in the sense 
that the surveys are sensitivity-limited, and it is likely that 
more evidence for mass loss could be found by increasing 
the sensitivity. Sonnett et al. (2011) hint at this possibil¬ 
ity by finding statistical evidence for comae in the averaged 
profiles of many asteroids, although not in any individual 
object. If real, their detection may point to electrostatic or 
gardening dust losses, or some other process acting broadly 
across the asteroid belt. 

Of the present sample, at least three active asteroids are 
likely to be driven by outgassing (133P, 238P, 313P), based 
on activity repeated in different orbits, corresponding to a 
fraction > 15%. The fraction of sublimation-driven objects 
is therefore of order / > 15 per million. The fraction of 
asteroids which contain ice, fice, could be much larger. We 
write fice = f / fd, where is the duty cycle, equal to the 
fraction of time for which the average asteroid loses mass. 
In the ice excavation model described earlier, fd is the ratio 
of the lifetime of an exposed ice patch to the interval be¬ 
tween exposures (caused by boulder impacts?) on the same 
body. The present level of ignorance about fd allows the 
possibility that fiee = 1, meaning that potentially all aster¬ 
oids could contain ice without violating the survey limits. 
This surprising result indicates that the dividing line be¬ 
tween comets and asteroids may be considerably less sharp 
than once supposed. 
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Table 1 

Summary of Orbital Properties 


Name 

T.f 

a'^ 

e® 


(7® 

O' 

(3200) Phaethon 

4.508 

1.271 

0.890 

22.17 

0.140 

2.402 

311P/PANSTARRS (P/2013 P5) 

3.662 

2.189 

0.115 

4.97 

1.936 

2.441 

P/2010 A2 (LINEAR) 

3.582 

2.291 

0.124 

5.26 

2.007 

2.575 

(1) Ceres 

3.309 

2.768 

0.076 

10.60 

2.556 

2.979 

(2201) Oljato 

3.299 

2.172 

0.713 

2.52 

0.623 

3.721 

P/2012 F5 (Gibbs) 

3.228 

3.004 

0.042 

9.73 

2.877 

3.129 

259P/Garradd (P/2008 Rl) 

3.216 

2.726 

0.342 

15.90 

1.794 

3.658 

(596) Scheila 

3.208 

2.928 

0.165 

14.66 

2.445 

3.411 

288P/(300163) 2006 VW 139 

3.203 

3.050 

0.200 

3.24 

2.441 

3.659 

(62412) 2000 SYi78 

3.197 

3.146 

0.090 

4.76 

2.864 

3.445 

P/2013 R3 (Catalina-PANSTARRS) 

3.185 

3.033 

0.273 

0.90 

2.204 

3.862 

133P/(7968) Elst-Pizarro 

3.184 

3.157 

0.165 

1.39 

2.636 

3.678 

176P/( 118401) LINEAR 

3.167 

3.196 

0.192 

0.24 

2.582 

3.810 

238P/Read (P/2005 Ul) 

3.152 

3.165 

0.253 

1.27 

2.364 

3.966 

P/2012 T1 (PANSTARRS) 

3.134 

3.154 

0.236 

11.06 

2.411 

3.897 

313P/Gibbs (P/2014 S4) 

3.132 

3.156 

0.242 

10.97 

2.391 

3.920 

P/2010 R2 (La Sagra) 

3.098 

3.099 

0.154 

21.39 

2.622 

3.576 

107P/(4015) Wilson-Harrington 

3.083 

2.638 

0.624 

2.79 

0.993 

4.284 


®^Tisserand parameter with respect to Jupiter 
'^Semimajor axis [AU] 

'^Orbital eccentricity 
'^Orbital inclination 
®Perihelion distance [AU] 

^Aphelion distance [AU] 
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Table 2 

Physical Properties 


Name 

lU 

Pv'^ 

pc 

5-1/® 

dm/dt ® 

(3200) Phaethon^ 

5-7 

0.08-0.17 

3.603 

0.58±0.01 

N/A 

311P/2013 P5 (PANSTARRS)^ 

<0.5 

0.29® 

7 

0.77±0.03 

7 

P/2010 A2 (LINEAR)® 

0.12 

0.1® 

7 

? 

N/A 

(1) Ceres® 

975 

0.090±0.003 

9.07 

? 

~6 

(2201) Oljato® 

1.8 

0.43±0.03 

7 

0.83±0.02 

5? (gas) 

P/2012 F5 (Gibbs)® 

1.8 

0.05® 

3.24±0.01 

- 

- 

259P/(Garradd)^ 

0.30±0.02 

0.05® 

7 

0.63±0.03 

<1.5 (gas), 0.01 

(596) Scheila® 

113±2 

0.038±0.004 

15.848 

0.71±0.03 

<3 (gas) 

288P/(300163) 2006 VW 139 ® 

3 

0.04® 

7 

? 

7 

(62412) 2000 SY178®° 

7.8±0.6 

0.065±0.010 

3.33±0.01 

0.64±0.03 

7 

P/2013 R3 (Catalina-PANSTARRS)“ 

<0.4 (multiple) 

0.05® 

7 

0.66±0.04 

7 

133P/(7968) Elst-Pizarro®® 

3.8±0.6 

0.05±0.02 

3.471±0.001 

0.65±0.03 

<0.04 (gas), 0.01,0.7-1.6 

176P/(118401) LINEAR®® 

4.0±0.4 

0.06±0.02 

22.23±0.01 

0.63±0.02 

0.1 

238P/Read®® 

0.8 

0.05® 

7 

0.63±0.05 

0.2 

P/2012 T1 (PANSTARRS)®® 

2.4 

0.05® 

- 

0.65±0.07 

- 

313P/2014S4 (Gibbs)®® 

1.0 

0.05® 

7 

0.72±0.02 

0.2 to 0.4 

P/2010 R2 (La Sagra)®'^ 

1.1 

0.04® 

7 

? 

4 

107P/(4015)Wilson-Harrington®® 

3.5±0.3 

0.06±0.01 

7.15 

7 

<150 (gas) 


'^Effective diameter (km) 

Geometric albedo 

‘^Rotation period 

‘^Color index. Solar color is B - V = 0.64±0.02 

^Inferred mass loss rate in kg s“^. Unless otherwise stated, the estimates are based on continuum measurements and refer to dust. N/A 
means that no mass loss rate can be specified because the loss is not in steady state, or for some other reason. 

^Value is assumed, not measured 

^Ansdell et al. 2014, ^Jewitt et al. 2013c, 2014c, Hainaut et al. 2014, ®Jewitt et al. 2010, 2011a^Kuppers et al., 2014, ®Tedesco et 
al. 2004, McFadden et al. 1993, Russell et al. 1984, ^Stevenson et al. 2012, Moreno et al. 2012, Novakovic et al. 2014, Drahus et al. 2015 
^Jewitt et al. 2009, ®Tedesco et al. 2002, Warner 2006, ®Hsieh et al. 201 Id, ^“Sheppard and Trujillo 2014, ^^Jewitt et al. 2014b ^^Hsieh 
et al. 2004, 2009a, 2011a, ^®Hsieh et al. 2011, Licandro et al. 2011, ^'‘Hsieh et al. 2011c, Moreno et al. 2013, Hsieh et al. 2013, 
O’Rourke et al. 2013 ^®Jewitt et al. 2015, Hui and Jewitt 2015, ^^Moreno et al. 2011b, Hsieh et al. 2012c, Hsieh 2014 ^®Veeder et 
al. 1984, Fernandez et al. 1997, Licandro et al. 2009, Urakawa et al. 2011, Ishiguro et al. 2011b. 
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Table 3 

Summary of Mechanisms’^ 


Name 

Sublimation 

Impact 

Electrostatics 

Rotation 

Thermal 

(3200) Phaethon 

X 

7 

7 

7 

/ 

311P/PANSTARRS (P/2013 P5) 

X 

X 

X 

/ 

X 

P/2010 A2 (LINEAR) 

X 

/ 

X 

/ 

X 

(1) Ceres 

/ 

X 

X 

X 

X 

(2201) Oljato 

7 

7 

7 

7 

X 

P/2012 F5 (Gibbs) 

X 

/ 

X 

/ 

X 

259P/Garradd (P/2008 Rl) 

7 

7 

7 

7 

X 

(596) Scheila 

X 

/ 

X 

X 

X 

288P/(300163) 2006 VW 139 

/ 

7 

7 

7 

X 

(62412) 2000 SYi78 

7 

/ 

7 

/ 

X 

P/2013 R3 (Catalina-PANSTARRS) 

7 

X 

X 

/ 

X 

133P/(7968) Elst-Pizarro 

/ 

X 

7 

/ 

X 

176P/( 118401) LINEAR 

/ 

7 

7 

X 

X 

238P/Read (P/2005 Ul) 

/ 

X 

X 

7 

X 

P/2012 T1 (PANSTARRS) 

/ 

X 

X 

7 

X 

313P/Gibbs (P/2014 S4) 

/ 

X 

X 

7 

X 

P/2010 R2 (La Sagra) 

/ 

7 

7 

7 

X 

107P/(4015) Wilson-Harrington 

7 

7 

7 

X 

X 


‘‘Z - evidence exists consistent with the process, x - evidence exists inconsistent with the process, ? - 
insufficient or mixed evidence exists 
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